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Abstract. We compute the magnetic susceptibility and specific heat of the spin-| Heisenberg model on 
the kagome lattice with high-temperature expansions and exact diagonalizations. We compare the results 
with the experimental data on ZnCu3(OH)6Cl2 obtained by Helton et al. [Phys. Rev. Lett. 98, 107204 
(2007)]. Down to ksT/J ~ 0.2, our calculations reproduce accurately the experimental susceptibility, 
with an exchange interaction J ~ 190 K and a contribution of 3.7% of weakly interacting impurity spins. 
The comparison between our calculations of the specific heat and the experiments indicate that the low- 
temperature entropy (below ~20 K) is smaller in ZnCu3(OH)6Cl2 than in the kagome Heisenberg model, 
a likely signature of other interactions in the system. 
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1 Introduction 

After many years of theoretical investigations, the nature 
of the ground-state of the spin- ^ Heisenberg model on the 
kagome lattice is still not known. Although all numerical 
studies have concluded to the absence of long-range mag- 
netic (Neel) order [1,2,3,4,5,6,7,8], many basic questions 
such as the existence of spontaneously broken symmetries, 
or the existence of a finite gap to magnetic excitations re- 
main open. In fact, many different states of matter have 
been proposed for the kagome Heisenberg antiferromag- 
net: Z2 gapped topological liquids [!),10], valence-bond 
crystals [11,12,1.3,14], critical spin liquids with gapless 
spinous [15, 12]. 

Recently, a promising spin- ^ antiferromagnetic insula- 
tor with an ideal kagome geometry, ZnCu3(OH)6Cl2, has 
been synthesized and studied for its magnetic properties 
[1 0, 17, 18, 19,20]. Because these studies did not detect any 
kind ordering (nor spin freezing) down to 50mK, it could 
represent one of the first and most remarkable realization 
of a 2D quantum spin liquid [21,22]. 

To extract some information about the low-energy physics 
of the kagome Heisenberg model from the experiments on 
ZnCu3(OH)6Cl2, it is important to first analyze in a quan- 
titative way the possible role of magnetic defects (and 
other "perturbations" to this model) in this compound. 
In this paper we compare the experimental data for the 
magnetic susceptibility x f^nd specific heat c„ obtained by 



Helton et al.[J^■] with calculations for the spin-^ Heisen- 
berg model on the kagome lattice based on exact diago- 
nalization (ED) data (partial spectrum of a 36-site cluster 
and full spectrum for 24-site and 18-sites clusters) and 
high-temperature series expansion [ '3,24]. Down to tem- 
perature ksT/J = 0.2, the experimental susceptibility 
Xcxp (T) can be very well fitted by that of the kagome lat- 
tice Heisenberg model with J ~ 190 K plus a contribution 
of about 4% of impurity spins with weak mutual inter- 
actions (modeled by a ferromagnetic Curie- Weiss temper- 
ature of ~ —6 K), likely mostly due to antisite disorder 
(Cu substituted to Zn on sites between the kagome planes) 
[-'0,25]. The low temperature specific heat is dominated 
by impurities (and other perturbations) below 2 K and 
by phonons above 15 K. In the intermediate range, the 
calculated specific heat appears to be larger than in the 
experiment. We comment on this feature at the end of the 
paper. 



2 Uniform static susceptibility 

The spin-i Heisenberg model reads: 



V.^ j'^Si - Sj - g^J^BH 



Sf 



(1) 
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where the sum runs over pairs of nearest neighbor sites on 
the kagome lattice and H is an external magnetic field. 
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To fix the notations, we define the (zero-field) uniform 
susceptibiHty per site x{T) as x{T) = ^ ^^^mp- 
where N is the total number of spins. 

The high-temperature expansion of x been com- 
puted up to order O (T^i^) (included) by Elst ner and 
Young [23]: 

x{T) = -fxty.it - ksT/J) 

Xth(i) = (l/4)t-i-(l/4)t-2 + --- (2) 

where Cq = 0.25(5/15)^ is the Curie constant and t the re- 
duced temperature. The truncated series Xn{t) = Yl^=o '^i^' 
at order n = 14 and n = 15 agree with a relative error 
smaller than 10^^ for t > \. Down to this temperature, 
they already provide good approximations to Xth(i)- The 
convergence of high-T series can be improved using Pade 
approximants (PA). In the present case the PA provide a 
reliable estimate of Xth(i) at least down to t '-^ 0.5.^ One 
representative PA (numerator of degree 8 and denomina- 
tor of degree 7) is displayed Fig. 1. 

On a small enough system, it is possible to obtain 
by ED the full spectrum (2-^ energy levels for N spins) . 
We have done so for two 18-site and 24-site kagome clus- 
ters (with periodic boundary conditions). Then thermo- 
dynamic quantities can be computed exactly as a func- 
tion of T. For bigger systems, where one can still compute 
some eigenstates by ED, one may use the approximate 
method described in reference [-''^J to compute thermody- 
namic quantities, footnote In this method, one constructs 
the density of states in each symmetry sector from the 
exact low and high energy states obtained from ED and 
approximating in between the unknown part of the spec- 
trum with a smooth density of states. This smooth part 
is constructed so that the first moments of the density of 
states (Tr[7J"]), in each symmetry sector of the finite clus- 
ter, are exact up to n = 5. This Ansatz guaranties that the 
thermodynamics becomes exact at low T (when thermal 
excitations only involve the eigenstates computed exactly) 
as well as at high T (when a re-summed high-temperature 
expansion up to is valid). 

The results for the susceptibility x(i) are shown in 
Fig. 1. Above t = 0.2, the relative difference between the 
(exact) TV = 18 and = 24 curves is smaller than 0.5%. 
We therefore make the (rather safe) assumption that our 
finite-size results are good approximations to the thermo- 
dynamic limit down to i„iin — 0.2. This represents a small 
gain over the coupled-cluster expansion of reference [27], 
which is valid above t ^ 0.3. 

The x(t) obtained with the approximate method for 
A^ = 36 sites also agrees (with a relative error smaller than 
2%) with the TV = 24 results down to t ~ 0.2. Shghtly be- 
low, the 36-site susceptibility is still increasing and might 
be a better approximation to the infinite-size limit than 
the 24-site curve. Still, it is not possible decide at which t 
finite-size (and/or errors due to the approximation in the 

^ Comparing the PA and the exact curves for 18 and 24-site 
clusters shows that the PA is in fact correct down to T 0.4J, 
see Fig. 1. 



density of states) will become too important. Safely, we 
only use the theoretical results (noted Xth) above i ~ 0.2 
to fit the experimental data Xcxp for the susceptibility. 

We fit Xcxp to a sum of contributions from the kagome 
spins Xs and the impurities Ximp in the following way: 

4C 

Xexp(r) = — Xs(fcsr/J) -I- Ximp(r) (3) 
with Ximp(T) = (4) 

where J is the (unknown) magnetic exchange in between 
the spins in the kagome planes, x an impurity concentra- 
tion, Co = 0.25{gfj,B)^ /ks the Curie constant and ^imp 
the Curie- Weiss temperature of the system of impurities. 
Eq. 4 is the leading term in a high-temperature expan- 
sion for the system of impurities, which provides a sim- 
plified picture of their interactions. To be applicable, T 
should therefore be large compared to 6'inip. We assume 
that the system of impurities does not perturb the the 
kagome spins. 

We optimized numerically the parameters so that Xs 
fits the theoretical results Xth for t > 0.2. As can be seen 
on Fig. 1 , an almost perfect agreement can be obtained be- 
tween Xs = 4^(Xexp(T)-Ximp(T)) (squares) and the the- 
oretical estimates for the kagome susceptibility xth with 
the following parameters: Co — 0.504 K cm'^/(mol of Cu) 
(equivalent to a gyromagnetic factor g — 2.32) , J — 190.4 K, 
X = 0.03655 and ^imp = —6.1 K. We note that the value 
of J is in rough agreement with the values reported in 
reference [IG] (17meV~200 K) and [27] (170 K). We also 
check a posteriori that the lowest temperature of the fit 
{t = 0.2 ~ 38 K) is much bigger than j^impj, so that a 
Curie- Weiss approximation for the impurities is justified. 
6 K is also approximatively the transition temperature re- 
ported in Zn2:Cu4_a;(OH)eCl3 for x < 0.6 (replacing some 
Zn by Cu between the kagome planes) [17], so this energy 
scale may correspond to some couplings for spins located 
between the planes, where magnetic impurities could sit. 
We eventually notice that down to i = 0.15 (that is below 
the lowest temperature used for the fit) , Xexp — Ximp con- 
tinues to increase and to follow the A^ = 36 curve. This 
suggests that the maximum of the kagome susceptibility 
could indeed be below t = 0.15. 

Rigol and Singh [,; . ] analyzed the same experimen- 
tal data with another high-temperature method (coupled 
cluster expansion) and obtained a somewhat different con- 
clusion. They argued that Dzyaloshinskii-Moriya (DM) in- 
teractions provide a better description of the low-temperature 
increase of the susceptibility than impurities. Although 
we agree that DM interactions are certainly present in 
ZnCu3(OH)6Cl2 and that they should affect the physics 
of the system (at least at low-temperatures), it is also 
clear that a few percent of impurities must be present too 
and should have a visible effect on the susceptibility, even 
at rather high temperatures. According to reference [27], 
free impurities cannot explain the sharp increase of Xoxp 
below 60 K. In our analysis, this issue is solved by al- 
lowing for a small ferromagnetic Curie- Weiss temperature 
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Fig. 1. (color online) Magnetic susceptibility per spin as a 
function of temperature. Dashed (magenta) curve : experimen- 
tal data (Helton et al.), multiplied by J/(4Co) as a function 
of ksT/J with J ~ 190 K and Co = 0.504 K cmV(mol of 
Cu). Black squares: Xs = 45^(Xexp(T) - Ximp(r)), obtained 
from the experimental data Xcxp by subtracting the contribu- 
tion Ximp of a concentration x = 0.03655 of impurity spins 
with a Curie- Weiss temperature ^imp = —6.1 K. Red (resp. 
cyan) curve : Exact xth for a A'^ = 24 (resp. 18) site kagome 
cluster with periodic boundary conditions. Blue curve: Results 
for A'^ = 36 spins obtained with the (approximate) method of 
reference [26]. Green curve: Fade approximant from the high- 
temperature expansion at order t~^'^ . The Fade approximant 
is not converged below t ~ 0.4 whereas the finite-size curves 
are practically converged to the thermodynamic limit down to 
t ~ 0.2. Below t = 0.2J, the later curves are only indicative. 

^inip — —6 K for the impurities.^ One can indeed see that, 
once Ximp has been subtracted, the experimental results 
show a saturation of x around 20 K, in rough agreement 
with the measurements of reference [l!)]. The location of 
the maximum we obtain is however quite sensitive to the 
value X of the impurity concentration. 



3 Specific heat 

The experimental data for the specific heat Cy{T) are only 
available at very low temperature where the size effects 
on the ED results are large. We therefore also applied 
the high-temperature entropy method [2H]. It combines 
three pieces of information about the system: 1) The high 

^ The Curie- Weiss temperature is an average of the differ- 
ent exchange constants. However, due to the complexity of the 
interactions between impurities (disorder), the ferromagnetic 
sign of Simp does not necessarily imply that they behave ferro- 
magnetically at low temperatures. 



temperature series expansion of Cy{T), up to T"^'' [23, 
24]. 2) The ground-state energy per site eg of the Hamil- 
tonian. Here we use the following estimate cq = 2{Q\Si ■ 
Sj\Q) = —0.44 [ ]. 3) The exponent a describing the low- 
temperature behavior of the specific heat: Cy{t ^ 0) ~ t". 
The method then provides a set of Cv(t) curves (for dif- 
ferent Pade approximants) which all satisfy exactly the 
following properties: i) c^{T ^ '-^ T°' , ii) Cy{T oo) ~ 
the series expansion, iii) Cy{T)dT = —Neo, and iv) 

Cy{T)/TdT = Nkg ln(2). When the value of eo and a 
are both correct , one usually gets a large number of very 
similar curves but if either is incorrect, only a few and scat- 
tered curves will be obtained (more details in Refs. [28, 
24]). In the case of the kagome antiferromagnet, the value 
of a is not known and the entropy method gives a reason- 
able convergence for a = 1 and a = 2 [24]. Motivated by 
the experimental observation [16] of a low-temperature Cy 
with an exponent smaller than one, we also include here 
a calculation with a — 0.5. 

The results for some valid'^ PA from orders 14 to 17 are 
displayed in Fig. 2 together with the exact result obtained 
by ED of the full spectrum of a 24-site kagome cluster. For 
t > 0.3, these results are practically exact. 

In the two scenarios a — 1 and 2 there is a significant 
dispersion of the results for t < 0.3 from one PA to another 
[:' !].'' In both case, c„ show a low-t peak or shoulder as 
found from ED of finite-size systems. The choice a = 0.5 
leads to a smoother Cv{t) and improves significantly the 
convergence. It is however not clear if this improved con- 
vergence for small values of a (0.5 a < 1) indicates that 
a is actually smaller than one or an artifact of the present 
entropy method, which might be "slower" to stabilize a 
low-t peak (as with a — 2, see Fig. 2), than a smooth 
curve (as with a = 0.5). In any case, this is clearly re- 
lated to the unusually large low-temperature entropy of 
the kagome system. 

Fig. 2 also displays the experimental results (black 
squares) obtained by Helton et al. [Ki]. The only parame- 
ter here is the exchange constant, taken to be J ~ 190 K 
from the fits of the susceptibility data. Above 15 K, the 
phonon contribution is dominant and we have to focus on 
the lower temperatures to analyze the magnetic contribu- 
tion. Below 15 K the order of magnitude agrees with our 
calculation but there is no quantitative agreement. Several 
"perturbations" such as weakly interacting magnetic im- 
purities or magnetic anisotropics should indeed contribute 
to the specific heat at such low temperatures. 

Results for the integrated entropy S{T) = Cy{x)/xdx 
are plotted in Fig. 3, where one sees that the choice of the 

By the entropy method, the entropy s(e) is obtained as the 
power a/ {a -|- 1) of a rational fraction (FA) of the energy per 
site e. The specific heat curve is then obtained parametrically 
through T(e) = l/s'(e) and c„(e) = - s' {ef / .s" [e) . Only the 
FA which satisfy s(e) > 0, s'(e) > and s"(e) < in the range 
]eo,0[ are physically "valid". 

* This is due to the finite order in the high-temperature 
expansion. Still, for a given value of a, we believe that this 
method gives a qualitatively correct picture for c„(r), even at 
low T. 
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Fig. 2. (color online) Specific heat c„ vs temperature T. 
Black squares: experimental data from Helton et al. [Ui], as- 
suming J ~ 190 K. Red triangles: exact specific heat of a 24- 
site kagome cluster. Green (full), blue (dashed) and magenta 
(dotted) curves: Ci, calculated by the entropy method for three 
different values of the low-temperature exponent (a = 2, 1 and 
0.5). All valid PA from order 14 to 17 with numerator and de- 
nominator of degrees > 3 are shown. For a given a and at each 
temperature, their dispersion provides a rough estimate of the 
error bars. 



low-T exponent a of Cy has practically no influence on the 
theoretical S{t = ksT/J) above t ~ 0.06 and that, around 
this temperature, the experimental value is significantly 
lower than in our calculations. Of course, subtracting the 
contribution of the phonons (hard to estimate quantita- 
tively) and from the impurities would make the discrep- 
ancy even larger. Concerning the impurities, one sees in 
Fig. 3 that an applied magnetic field of 5 and 14 Teslas 
is enough to significantly reduce S{t) for t < 0.06. Such 
fields are low in comparison to J but of the order of the 
estimated coupling between the impurities. The difference 
between the curves at and 5 (or 14) Teslas may thus 
provide a rough estimate of the contribution of the impu- 
rities. Note that S{t) at 5 Teslas become closer to the S{t) 
computed by the entropy method with a ~ 2. The experi- 
mental value a < 1 could be due to the impurities and the 
actual value of a for the kagome spins could be larger than 
1. In any case, at < = 0.06(~ 11 K), the experimental en- 
tropy (~ 0.2 In 2) is thus at least 7% of In 2 below the the- 
oretical estimates (0.27 In 2). This seems a rather robust 
indication that some additional interactions play some role 
in this energy range, by freezing some degrees of freedom 
of the spins in the kagome planes and pushing the corre- 
sponding entropy to higher energies. We may mention in 
particular DM interactions [27], non-magnetic impurities 
in the kagome planes ("dilution") [20] and interactions 
between impurities and the kagome spins. We conclude 
that 15^20 K is a minimal temperature for a kagome lat- 
tice Heisenberg model description of ZnCu3(OH)gCl2 to 
be valid. 
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Fig. 3. (color online) Entropy 5* vs temperature. Square, tri- 
angles and circles: experimental data from Helton et al. [Ki] 
with magnetic field B= 0, 5 and 14 Teslas, plotted as a func- 
tion of fcsT/J with J = 190 K. Green (full), blue (dashed) 
and magenta (dotted) curves: entropy calculated by the en- 
tropy method for a = 2, a = 1 and q = 0.5 (same PA as in 
Fig. 2). 
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